Entanglement of two blocks of spins in the critical Ising model 



o 
o 

(N 



(N 



Oh: 



> 
O 

o 
p 

00 

o 

oo 
O 



X 



p. Facchi,^'^ G. Florio,'^'^ C. Invcrnizzi,^ and S. Pascazio'^'^ 

' Dipartimento di Matematica, Universita di Bari, 1-70125 Bari, Italy 
^INFN, Sezione di Bari, 1-70126 Bari, Italy 
Dipartimento di Fisica, Universita di Bari, 1-70126 Ban, Italy 
Dipartimento di Fisica, Universita di Milano, 1-20133 Milano, Italy 

We compute the entropy of entanglement of two blocks of L spins at a distance d in the ground 
state of an Ising chain in an external transverse magnetic field. We numerically study the von 
Neumann entropy for different values of the transverse field. At the critical point we obtain analytical 
results for blocks of size L = 1 and L = 2. In the general case, the critical entropy is shown to be 
additive when d ^ oo. Finally, based on simple arguments, we derive an expression for the entropy 
at the critical point as a function of both L and d. This formula is in excellent agreement with 
numerical results. 
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I. INTRODUCTION 

A comprehension of the features of entanglement in 
systems with many degrees of freedom, such as quantum 
spin chains, is currently one of the most challenging prob- 
lems, at the borderline of quantum information science 
[H and statistical physics. In the last few years several 
measures of entanglement have been proposed 0, H, Q 
and calculated (analytically in the simplest cases, oth- 
erwise numerically) for the ground states of many-body 
systems @. 

Despite accurate investigations and different propos- 
als, there is still no consensus on the correct character- 
ization of the multipartite entanglement of the ground 
state of a many-body system. We will consider here the 
entanglement entropy, a measure that can sometimes be 
tackled by analytic investigations and for which quantum 
field theoretical methods can be employed. The entangle- 
ment entropy is just the von Neumann entropy associated 
with the reduced density matrix, that is the entropy of a 
subsystem of the chain, and was explicitly evaluated for 
quantum spin chains [j, 0, H, lol. [Tol. fill . 

One of the most striking features of the entanglement 
entropy is its universal behavior at and close to a quan- 
tum phase transition. Indeed, it is found that entropy in 
non-critical systems generally tends to saturate towards 
a finite value as the size of the subsystem increases, but 
this value (logarithmically) diverges with the size of the 
subsystem as the system approaches a quantum critical 
point. Close to a quantum critical point, where the cor- 
relation length <^ is much larger than the lattice spacing, 
correlations are described by a 1-1-1 dimensional quantum 
field theory and at the critical point, where ^ diverges, 
the field theory is also a conformal field theory [l2]. In 
the latter case, the behavior of entropy calculated by ana- 
lytical and numerical techniques for several spin systems 
is confirmed by the predictions of the corresponding field 
theory. 

In this work we extend the characterization of the 
entanglement entropy to a more general subsystem, in 
which the correlations between two disjoint blocks of 



spins and the rest of the chain are studied as a function 
of the distance d between the blocks and their common 
size L. This entropy of entanglement will be denoted 
S(L, d) and will be analyzed by analytical and numerical 
methods. 

The physical system we shall consider is the Ising 
model in a transverse magnetic field, since it fulfills a 
convenient combination of requirements. It is solvable, 
its ground state can be computed by using well-known 
analytical and numerical techniques [7[ and, at the same 
time, it successfully describes a rich spectrum of physi- 
cal phenomena, that include the ordered and disordered 
magnetic phases, connected by a quantum phase transi- 
tion [ll|. 

We will analytically compute the entropy of entangle- 
ment at the quantum critical point (QCP) for blocks of 
L = I and L ~ 2 spins, and will tackle the problem nu- 
merically for larger values of L. We will first study the 
behavior of the entropy as a function of the magnetic 
field A, then at the QCP, A = 1, as a function of the 
distance d between the blocks and their size L. We will 
investigate the limits d — > and d — > oo. Our results 
will include as a particular case (d = 0) the logarith- 
mic behaviour of the entropy of a single block of L spins 
at criticalit y S l = h log L -\- /C, where /C is a constant 
0, 0, i, i,M im, Q [il, m. We will also show the 
additivity of S{L, d) at the critical point as d —f oo. Fi- 
nally, we will plot S{L, d) as a function of both L and d, 
getting an accurate idea of the features of the entropy at 
the QCP. 

This paper is divided in six sections. In Section 
n we review previous works on spin chains, following 
IJ, 0, [13, E H, m , where the ground state of the Ising 
model is computed: the explicit expressions obtained will 
be used in the following sections in order to obtain the 
reduced density matrix and the Von Neumann en- 
tropy Sl of L contiguous spins. In Section [llll we extend 
the definition of the correlation matrix given in 0, 0| 
to a bipartition of two blocks of L spins separated by a 
generic distance d. We define here the entropy S{L,d), 
describing the entanglement of the two blocks with the 
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rest of the chain. By making use of the newly defined re- 
duced density matrix QL,d^ in Section [TVl we analytically 
compute the entropy for blocks of one and two spins at 
the critical point. In Section |V] we carry out numerical 
computations of S{L, d) for several sizes of the blocks. 
Finally, we plot the entropy as a function of the size L 
of the blocks and their reciprocal distance d, in order 
to get a general idea of the features of the entropy of 
entanglement at the critical point. Our results are sum- 
marized and discussed in Section IVII In the Appendices 
we included, for self consistency, additional material and 
explicit calculations. 



II. GROUND STATE OF THE ISING MODEL 

The Ising chain in a tranverse field consists of 2N + 1 
spins with nearest neighbor interactions and an external 
magnetic field, described by the Hamiltonian 



Hi 



-J 



E 

-N<i<N 



(1) 



Here i labels the spins (we take an odd number of spins 
for simplicity), J > and we consider open boundary 
conditions, cr^_|_i = 0. erf (/i = x,y,z) are the Pauli 
matrices acting on spin i. The determination of the 
ground state proceeds with the Jordan- Wigner transfor- 
mation in terms of Dirac or Majorana fermionic opera- 
tors [13, m, E, [iSl • Here it will be convenient to consider 
Majorana fcrmions, whose operators are defined by 



a2/-i 



n 

\m<l 



a2l 



^r, (2) 



with —N < I < N. They are hermitian and obey anti- 
commutation relations, 



Figure 1; (Color online) Two blocks of L adjacent spins at a 
distance d. The state gL,d is obtained from the ground state 
\tf)o) of the spin chain by tracing out the spins that do not 
belong to the blocks. 



such an invariance: {d2mO'2n} = {d,2m-i6.2n-i) = 0,'im,n 
with m 7^ n, while (a2m-i^2n) = idm-n depend only on 
the difference m — n. Therefore, the block correlation 
matrix becomes independent of k and reads 



/ Ho 

Hi 



n_i 

Ho 



L + l 



\ 



Hn 



(7) 



with 
Til = 

(8) 

where the real coefficients gi are given, for an infinite 
chain, by 



-i(a2i-iao) 






i{a2id-i) 




) 



1 

Ztt 



(cos (/) — X) + i sin cj) 
(cos (j) — A)2 + sin^ (f> 



(9) 



Thus, Tf^ is a real, skew-symmetic 2L x 2L matrix, i.e. 
(^l)^ = — F;^, since all the blocks H; have the property 

(H,)^ = -n_,. 



— 1 {'^mj^n} — '^^mn^ (3) 

and their expectation values in the ground state {ipo) 

(l/'O I Oman I l/'o) = (dmCLn) = 6,nn + i^^mi (4) 

with —2N — 1 < ni,n < 2N, completely characterize 
iV'o)- Consider now a block of L contiguos spins labeled 
by i with 



k<i<k + L-l, 



(5) 



with k > —N and k + L — 1 < N. The expectation values 
of the Majorana operators of the block are encoded in the 
2L X 2L submatrix 



^ ( (^rn ) ^rnn ) 5 



(6) 



with 2k - 1 < m,n < 2k + 2L - 2. 

We arc interested in the thermodynamic limit of an 
infinite chain, TV —f oo. In such a limit, the ground state 
becomes translation invariant, and all correlations inherit 



III. ENTROPY OF TWO BLOCKS OF SPINS 

The entropy of a single block of L contiguous spins in 
the critical regime can be obtained by very accurate nu- 
merical results H, 01 and analytical conformal field theory 
calculations Given the ground state lipo), one finds 

the reduced density matrix 



QL = Tr^L{|'0o)(V'o|}, 



(10) 



where the trace is over all spins that do not belong to the 
block, and its entropy 



Sl = -TrlgLloggL}. 



(11) 



For definiteness, in this paper we will fix the base of log- 
arithms to 2. The calculations of Refs. [1, 0|, yielding an 
expression of the reduced density matrix qi^ and the en- 
tanglement entropy Sl of L adjacent spins in the ground 
state iV'o) are reviewed in Appendix |^ The key point 



3 



is the following: the block entropy (fTTj) is the sum of L 
terms 



1=1 



l + vi 



where 



H{x) ~ — .xlogx — {1 — x) log(l — x) 



(12) 



(13) 



is the Shannon entropy of a bit, and ±w;, with \ < I < 
L, are the pairs of (purely imaginary) eigenvalues of the 
block correlation matrix of Eq. Q . 
In the continuous limit 



Sl = - log L + JC{L), 
6 



(14) 



where 1/6 derives from the central charge c = 1/2 of a 
free massless fermionic field and 



/C(L) ^/C + O ( ^ 



(15) 



with L<x<L + d in the (2L + d) x {2L + d) matrix 



2L+d 



2L+d 



A 



A 



id) 



A 
A 



-L-d 
{d,L) 



L+d ^ 



(L.d) 



A 



(19) 



For example, let us consider the case of two blocks of 
L = 2 spins at a distance d = 3. In this case the 7x7 



^ 7 



A 

2L+d ^ 




reads 
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(20) 



JC being a constant. In this section we want to extend 
this approach and construct the density matrix gi^ ^i of a 
subsystem of two blocks of L adjacent spins situated at 
a distance d, studying their entanglement with the rest 
of the chain. Sec Fig. [TJ To this aim, one starts by com- 
puting the matrix of a single block of adjacent spins 
([7]), and then traces out the central d spins as follows. 

We define the 4i x 4L correlation matrix F;^ ^ of two 
blocks, each of L spins, situated at a distance d {d, like L, 
are expressed in units of the distance between adjacent 
spins, and arc therefore dimcnsionless) 



^ L,d 



' ^0 ^-L-d 



\ ^L+d ^0 



A 



-L-d 



(16) 



where Ai^'^ = ^x"^^ with 



n^+l Ha; 



na-_A/+2 



(17) 



The matrix ASt^ = ^— A^^^ is a Toeplitz matrix, and 
^ has the property 



pA 



2L; 



(18) 



i.e., when the distance between the two blocks is zero, 
F;^ ^ becomes equal to the matrix ([7]) of a single block 
of size 2L. The matrix ^ in Eq. is obtained by 
tracing out the d rows and d columns that are labeled 



and we have to cancel the columns whose first element is 
labelled by —2, —3, —4 and the rows whose first element 
is labelled by 2, 3, 4, obtaining 



pA _ 
^ 2,3 — 



Ho n_i n_5 n_6 

Hi Ho n_4 n_5 

Hs n4 Ho n_i 

He Hs Hi Ho 



(21) 



that is again a real skew-symmetrix matrix. 

The entanglement of the two blocks of spins reads 



S[L,d) = -TiigL.dlogQL.d), 



(22) 



where g^^d is the density matrix of two blocks of L adja- 
cent spins at a distance d. Exactly as for a single block, 
the above entropy can be given an explicit expression in 
terms of the eigenvalues ±wj, with 1 < / < 2L, of F^^, 
analogous to 



2L , 

S{L,d) = Y^H[- 



1=1 



n 



(23) 



with H given by 

Before investigating the behavior of Eq. ((22|) , it is in- 
structive to look first at some simple examples. 



IV. ANALYTICAL RESULTS 

A. Entanglement of Two Single Spins 

We consider the Ising chain in a critical transverse 
magnetic field Ac = 1. At A = Ac, the coefficients gi of 
the reduced correlation matrix F^ ^, defined in Section 



4 



mil and given by Eq. can be computed analytically, 
yielding 



91 



2tt 



27r 



7r(/ + i)- 



(24) 



In order to compute the entanglement entropy S{L, d) 
of two single spins at a distance d, we have to calculate 
the eigenvalues of the correlation matrix ^, where L = 
1 and d> 0. Equation (fT6|) reads 



0{l/d'^) coincide in the limit d — !■ +oo, we have 



- log — + log 



More on this phenomenon later. 



(30) 
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-2 

7t(2;+i) 




ir(2/+l 



2 



7r(-2i+l) \ 


-2 







,(25) 



where I = d+1. The eigenvalues are the solutions to the 
characteristic equation 



det(rf , - /i) = 0. 



(26) 



For L — 1, there are four eigenvalues zLivi and ±W2, 
with 



2 - 1)2 +4^2 ± 1 

TT 4P - 1 



(27) 



The eigenvalues of the reduced density matrix gL^d are 
(1 ± !^fc)/2 and the von Neumann entropy reads 



S{l,d)^H 



H 



1 + 1^2 



(28) 



where H is given by (|T3l) . 

It is interesting to consider the cases d = and d ^ 
1. In the former case we have ^1,2 = 2(-\/T3± l)/(37r), 
whence 



5'(1,0) = F 



13 + 2 



H 



13 



37r 



1, f I 16 - 77r2 

-2^°ni^+^^ 



log 1 



137r 



1 



16-4x/T37r + 37r2^ 
Stt 



■ log IH - 

37r ^ V 37r2 - 47r - 16 



(29) 



In the latter case, since the eigenvalues 1^1.2 — 2/tt 



B. Entanglement of two blocks of L — 2 spins 

If L = 2, we have to compute the eigenvalues of 
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Hi 
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-(i+d) 

Ho 
Hi 
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Ho 
Hi 



-(2+0 

(1+0 



n-1 
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(31) 



with I = d+1. The characteristic equation ((26l) is of 
8th degree, but can be reduced to a quartic equation in 



t 



r + pr + qr + rt + s = o, 



(32) 



which has an exact solution. The coefficients p, q, r and s 
are functions of the distance d and are explicitly written 
in Appendix [Bl 

The matrix F^ ^ will have the eight eigenvalues ±/ifc = 
ztivk, with k = 1,2,3,4, two for each (negative) root of 
t. The eigenvalues of the reduced density matrix j32,d are 
(1 ± i/fc)/2 and the entropy reads 



S{L,d) = Y.^ 



k=l 



(33) 



with H given by (fO|) . 

In Fig. [2] we plot the eigenvalues i^k versus the dis- 
tance d. Note that the eigenvalues quickly saturate at a 
distance d ~ 5 between the blocks. This means that the 
entanglement between the two L ~ 2 blocks reaches its 
asymptotic value for rf > 5. 

The asymptotic values of the eigenvalues are solutions 
to the equation obtained by taking the limit d — > 00 of 
([5^ . One gets 



f . f . 2 \/l3 + l 

t^l 00 = 1^3 00 = , 

TT 3 



1^2(00) = (^4(00) 



2 ^/I3- 1 



(34) 



in agreement with Fig. [21 Note that they coincide with 
(P7)) evaluated at I = 1, i.e. with the eigenvalues of 



5 



n±i 



(cti T «a-2)/2, Eqs. (HH) and (HIl) greatly 
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Figure 2: (Color online) Eigenvalues of the reduced density 
matrix of two blocks of L = 2 spins versus their distance 
d at the critical point \c = 1. Note that the v^s reach a 
saturation value at d ~ 5. Here and in the following figures, 
d is expressed in units of the distance between adjacent spins, 
and is therefore dimensionless. 



the reduced density matrix of two spins at a distance 
d = / — 1 = 0. Therefore, for d — > oo the eigenvalues 
coalesce into pairs and the spectrum of two blocks (of 
two spins) coincides with the spectrum (with degeneracy 
2) of a single block. This phenomenon, which implies 
the asymptotic additivity of block entropy, is indepen- 
dent of the blocks dimension L and will be discussed in 
full generality in Sec. IVB2I 



V. ENTROPY OF TWO BLOCKS OF SPINS 

In this section we look at the entanglement entropy 
S{L,d) of two blocks of spins when the magnetic field 
varies. At the critical point, Ac — 1, we will find an ex- 
pression for S{L, d) in terms of the entropy Sl of a single 
block, given in Eqs. (|lip - (fT4|) . and investigate its limits 
d —I- 0, d oo. We will combine numerical estimates 
with analytical methods. 



A. Entropy versus A 

We start by evaluating the entanglement entropy 
S{L, d) versus the magnetic field A. The entanglement 
between two single spins at a distance d and the remain- 
ing part of the chain was investigated in Refs. [23l. [2^. [25j 
as a function of the magnetic field A. A generalization to 
a comb of m spins, spaced d sites apart can be found in 
[26j . We now generalize these results to the case of two 
arbitrary blocks of L spins at a distance d. 

In general, the presence of a gap between the blocks 
yields a larger entropy for all values of the magnetic 
field A. Let us start examining the situation at zero 
magnetic field. At A = 0, from ^ one gets that 
gi = Si^-i. Therefore, since 11; = for / 7^ ±1 and 



simplify. In particular, for any L > 1, the off-diagonal 
blocks in (fTB]) read, for d = 0, 



/O 



A 



\0 .. 

and yield a tridiagonal block matrix 




(35) 



^ L.O — ^ 2L 



/ n_i 

n+i n_i 

y 







0/ 



while, for any d > 0, = - 

gets 



(36) 



and one 



(37) 



Thus, at d = 0, the characteristic polynomial is 

det(rf^o - ^) = dct(r^^ + if^~\ (38) 

hence vi =Q and i^; = 1 for 2 < Z < 2L, so that 

2L 

2 1'^"''' " V 2, 



1=2 



(39) 



On the other hand, for d> 1, one gets 

dot{rl, -f,)= det(rf - = /(/i^ + i)2L-2^ (40) 
hence = 1^2 ~ and vi ~ I for 3 < I < 2L and 

T 



SiL,d) = 2S{L,0) = 2H 



(41) 



This is intuitively clear: at zero transverse field S{L, d) 
follows exactly an area law [27j and in 1 dimension the 
presence of a gap doubles the area of the boundary, dou- 
bling the entropy. See Fig. [31 For nonzero values of the 
magnetic field, there are correction to the area law, due 
to correlations between the two blocks. An entropy in- 
crease is still natural, but it turns out to be smaller than 
the factor 2 that one would naively expect for a doubled 
boundary. We will show that the factor 2 can be recov- 
ered in the limit of large gap d: this is the phenomenon 
of asymptotic additivity of entropy mentioned at the end 
of Sec. HVBl See Sec. IVB] 

We show the results of some numerical investigations 
in Fig. [3] Let us first look at the case of small L (~ 2-;-3). 
For d > entropy is not maximum at the critical point 
A = 1, but rather for some A < 1. On the other hand, 
for bigger L, the maximum entropy is always reached at 
the critical point A = 1 and its value grows with the size 
L of the blocks. 



6 



One also notices that, far from criticality, the entropy 
has a very weak dependence on the value of the gap d > 
1, and thus an area law is a very good approximation. 
The largest deviations are at the critical point, when the 
correlation length diverges. At fixed L, the ordinate of 
the cusp at A = 1 is always an increasing function of 
d. We now turn to the study of the critical case and 
endeavor to find some interesting analytical expressions. 



B. Critical chain 

1. Blocks of contiguous spins (d —> Q) 

In the limit d ^ 0, the entropy of entanglement ((22] 
must reproduce the single-block result (jlip as a particu- 
lar case: 



S{L, 0) = S; 



2L- 



(42) 



This is a simple consistency check and was numerically 
verified when the magnetic field is critical, Ac = 1. As 
a byproduct, this enables us to obtain the value of the 
constant K, via the logarithmic fit 



S'(L,0) = ilog2L + /C 
6 



O 



We obtain 



/C = 0.690413, 



(43) 



(44) 



with an error ~ 9 • I0~^, corroborating the results in [1]. 
An accurate fit enables us to give a precise estimate of 
the corrections in l/£, but more on this later. 



2. Asymptotic additivity of entropy (d oo) 

The plot in Fig. 2] shows that in the limit d ^ oo the 
two-block entropy is accurately fitted by 

S{L,oo) = \\ogL + 2ICiL) ^2Sl. (45) 

Entropy becomes therefore additive at large distances d. 
The physical meaning of this result is that the entangle- 
ment entropy of two separated blocks of L spins becomes 
twice the entropy of a single block L when the distance 
between the blocks approaches oo, i.e. it becomes much 
larger than the size of a block. Therefore, at the crit- 
ical point, the quantum correlations between two finite 
blocks of spins saturate at a certain distance. We now 
explain this result, which turns out to be valid for every 
value of the magnetic field A. 
Note that from ^ one gets 



hm gi=Q, 



(46) 
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Figure 3: (Color online) Entropy of the reduced density ma- 
trix for blocks of (a) L = 2, (b) L = 3, (c) L = 8 and (d) 
L — 15 spins versus the external magnetic field A for several 
distances d between the blocks. Circles d = 0, squares d = 10, 
diamonds d = 50. Squares and diamonds are indistinguish- 
able in (a) and (b), and are barely distinguishable in (c) and 
(d) only around the critical value A = 1. 
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by Ricmann-Lcsbegue lemma. Thus, from 
lim Ui = 0, 



and from pT)) 



lim Ai^^ 

X — '■OO 



0, 



VL. 



(47) 



(48) 



A^x"^ accounts for the residual correlation of two blocks 
of spin at a distance x. 

Therefore, the limit of the reduced correlation matrix 
of two blocks ITCl) reads 



yA 

^ L.i 



a — * + C30 



lim 



A 



A 



A, 



{L) 



(L) 

-L-d 



(49) 



for all L. The two blocks become independent and the 
limiting spectrum of ^ is given by the spectrum of a 

single block Vj^, with degeneracy 2. As a consequence, 
the entropy (|23p becomes additive in the limit 



2L 



S{L, oo) = 



1=1 \ 



2Sl, (50) 



where ±Wi(r^^) with 1 < I < 2L denote the eigenval- 
ues of Tf^ ^ and ±Wm(r;^) with 1 < m < L denote the 

eigenvalues of T^. We stress again that these results are 
valid for all values of the magnetic field A. See the intro- 
ductory comments in Sec. IV Al The area law discussed 
there is restored for sufficiently large d when the correla- 
tions between the two blocks arc negligible and the two 
boundaries become "independent" . 



3. General behavior of S{L, d) at the critical point 

We now turn to the problem of describing the entan- 
glement of two blocks of L spin at an arbitrary distance 
d with the remaining part of the critical Ising chain. In 
order to find a function of L and d we make two assump- 
tions: we require that the entropy be a function of all the 
scales of the problem; moreover, the dependence must be 
logarithmic. 

Let therefore 

S[L,d) = i(21og(L-a)-21og(i-Fd) 

+ log(2L + d-a) + \og{d H- a) -f , (51) 




100 



Figure 4: (Color online) Saturation value of the critical block 
entropy (A = 1). The numerical values of S{L,<x) are fitted 
by ilogL-H2/C(L). 



with < a < 1 and € M. See Fig. [H The quan- 
tity a fixes the position of the end of each block, a = 
corresponding to the central point between two adjacent 
spins, while a = 1 to the position of the last (or first) 
spin. Clearly, a detects granularity in the chain and the 
corrections due to a will be important for small values of 
L and/or d. We obtain 

SiL,0) = ilog(2i-a)-Hilog(l-|) 

+ i(loga + /3) 
6 



ilog(2i)-f i(loga + /3), L 



whence 

On the other hand, 
S{L, +oo) 

whence 

In conclusion, 
1 



loga + /3 = 6/C. 



(52) 
(53) 



21og(L-a)+/3 
1 



-log(L)-f -/3, 



(3 = 12/C. 



(54) 



(55) 



S{L,d) = -[2\ogiL-a)-2logiL + d) 

+ log(2L + d~a)+ log{d + a) ~ 2 log , 

(56) 



with 



a = 2-^'^ = 0.0566226. 



(57) 
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Notice that there are no free parameters. Moreover, in 
the realm of vahdity of CFT, when d, L ^ 1, one gets 

S{L,d) ^ i(21ogi-21og(i + d) 

+ log(2L + rf) +logrf) +JC, 



(58) 

that agrees with the results of Calabrese and Cardy Q 
when one adds a missing addendum in their formula 
(3.32). 

From we get 



S{L,0)^-\og{2L)+JC{L), 



(59) 



where 

ICiL) = -6loga + -log(l--)+-log(l-- 



GO. 



(60) 



5 a 3 17 
12L " 16l2 ~ 144l3' 
This formula is in excellent agreement with numerical 
results. It provides the explicit expression in Eq. (I15p 
and was used in Eqs. ([13])- (03]) and Fig. [H See also the 
discussion at the end of Sec. IV B II 
The global behavior of 



S{L, d) 



i(21og(L-a)-21og(L + d) 
+ log(2L + d-a)+ \og{d + a) ^ 
i(21ogL-21og(L + d) 
+ log{2L + d)+logd) +IC{L), (61) 



with a and IC{L) given by ([57]) and (pO)) . respectively, is 
displayed in Fig. [S] The fit is accurate up to one part in 
10^ for small L (< 10) and one part in lO*" for L > 10. 
Notice the logarithmic L dependence for d = and the 
saturation effect for L/d <Si 1. A section of Fig. O is 
displayed in Fig. [6] In particular, the inset shows the 
asymptotic behaviour of the entropy and its saturation. 




Figure 5: (Color online) Critical entropy (A = 1) between two 
blocks of L spins at a distance d and the remaning part of the 
(infinite) chain. 
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Figure 6: (Color online) Critical entropy (A = 1) between 
two blocks of 1/ = 30 spins and the remaning part of the 
chain versus d. Inset: L = 30 with < d < 140. 



4. Behavior of the critical entropy for small d 

Both in Figs. [5] and [6] one notices for all values of L 
a sharp entropy increase at small values of d from d = 
to d = 1. This corresponds to the two ordinates of 
the cusps in Fig. [31 Let us endeavor to interpret this 
phenomenon on the basis of the formulas derived in this 
section. Equation ([55)1 yields 



AS* 



S{L,l)-S{L,Q) 



1 
6 



-2 log 



L + 1 



log 



2L~a 



log ■ 



, 1 + a 

log 

6 a 



log a = /C 



(62) 



where we used Eq. (1571) in the last equalities. This agrees 
very well with Fig. |6l and explains why AS* is largely 
independent of L in Fig.[5l More to this, the final result 
in Eq. (|62p yields a suggestive interpretation of the fitting 
parameter a in Eq. (|57p and of the constant /C: they turn 
out to be related to the entropy increase AS* associated 
with the opening of a d = 1 gap (one qubit) in an interval 
of 2L contiguous spins. It is therefore not surprising that 
K, = —(1/6) logo:, being 0{AS), be also necessarily of 
order 1. 



VI. CONCLUSIONS 

We provided an analytic and numerical treatment of 
the entanglement entropy of two disjoint blocks of spins 
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as a function of their length and distance in the quan- 
tum Ising model with a transverse magnetic field. We 
gave an analytic expression of the entropy at the criti- 
cal point, for two blocks of 1 and 2 spins at a generic 
distance. We showed that the presence of a gap always 
yields an entropy increase. At criticality, due to a log- 
arithmic correction to the area law, this increase is less 
than a factor 2 for all values of d, and becomes 2 for 
large d, when asymptotic additivity takes place. We also 
showed that, interestingly, the entropy of the two blocks 
can be written in terms of the entropy of a single block of 
spins, that, at the quantum phase transition, grows loga- 
rithmically with the size of the block. We have also given 
an accurate idea of the general features of the entropy of 
two blocks as a function of their size and distance. 

The behavior of the entanglement in a critical spin 
chain agrees with well known results in conformal field 
theory, where the geometric entropy (analogous to the 
spin block entropy, but defined in the continuum) can be 
computed for 1 + 1 dimensional theories [13, [iBl- The 
translation of field theoretical methods and ideas in the 
language of quantum information will hopefully enable us 
to make use of additional results for an arbitrary number 
of disjoint intervals [1, [2l|. This would of great inter- 
est from the point of view of multipartite entanglement. 
The study of the statistical distribution of bipartite en- 
tanglement for different bipartitions [2^ is a useful tool 
for the analysis of multipartite entanglement. A deeper 
comprehension of the dependence of entropy on distances 
and sizes of blocks of spins could yield information about 
the role of quantum phase transitions in the generation 
of multipartite entangled states [1^, [13] . 

Acknowledgments 



Then, V defines a new set of Majorana operators, c,]„ = 

Cm 5 

2L 

Cm = ^ ^ ^m.n^ii-j (-^^) 
n=l 

that satisfy the same anticommutation relations as the 
d„'s, and have correlation matrix r£. The structure of 
r£ implies that mode C2i~i is only correlated to mode 
C2/. In the language of fermionic operators, one gets L 
spinlcss fermionic modes 

. _ C21-1 + iC2l 

ci = ^ , 

{q. Cm} = 0, {cj. Cm} = Smn, (A3) 

that, by construction, fulfill 

(QCm) = 0, (c|cm) = <5;m^-^- (A4) 

Thus, the L (nonlocal) fermionic modes are uncorrelated, 
so that the reduced density matrix can be written as a 
product 

ql = Pi® ■ ■ ■® PL- (A5) 

Now, the density matrix pi, 1 < I < L, has eigenvalues 
(1 ± i'i)/2 and entanglement entropy 

Sipi) = H (y^P-) , (A6) 



We thank Matteo Paris for interesting discussions. 
This work is partly supported by the European Com- 
munity through the Integrated Project EuroSQIP. 



where H is the Shannon entropy of a bit (fT^]) . Therefore 
the spectrum oi gi^ results from the product of the spectra 
of the density matrices pi, and the entropy of is the 
sum of the entropies of the L uncorrelated modes. 



Appendix A: COMPUTATION OF Sl 

We review here the computation of the entropy of the 
reduced density matrix ql = Tr^i|'0o)('0o| for L adja- 
cent spins [1, [^l- In the limit of infinite chain (iV oo), 
a given finite section of the chain is fully translational 
invariant and describes the state of any block of L 
contiguous spins. The density matrix can be recon- 
structed from the restricted 2L x 2L correlation matrix 
of Eq. ([7]). In particular, a direct way to compute 
the spectrum of g^ and its entropy Sl from is the 
following. 

The matrix can be put into a block-diagonal form 
by an orthogonal transformation and its eigenvalues are 
purely imaginary and come in pairs, iii/j, and {i^il < 1, 
with 1< I < L. Let F G SO{2L) be the special orthog- 

is block-diagonal 



onal matrix such that r£ 



e 

1=1 



1 

-1 



(Al) 



1=1 



(A7) 



This is Eq. (|12p of the text. Summarizing, for arbitrary 
values of the magnetic field A and in the thermodynamic 
limit, A'^ — > oo, the block entropy Sl of the ground state 
of the Ising model is given by the sum (|A7p . where ±ii/j 
are the pairs of imaginary eigenvalues of the block corre- 
lation matrix Ff of Eq. © . 



Appendix B: EIGENVALUES OF F^^ 



The coefficients in Eq. (|32)) are functions of the dis- 
tance d between the blocks and read 
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Pid) 

q{d) 

rid) 
sid) 



2^(3 + 2d)4(5 + 2d)^ 

212(3 + 2d)^(5 + 2d)^ 
7r4A 



5303 



24314B 41528^2 101445^ 



203297 + 391466B- 



9 9 
2841841^2 
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3652160^3 2617216B* 3299845^ 171525^ 



2^2(2 + 

TT^A 
232 



12271 



27 81 81 81 

68116S 158795^2 198074^^ 1390005* 173125^ 



3 



27 



81 



81 



(l + d)*(2 + d)«(3 + d)^ 



8965^ 
81 



(Bl) 



where A = {I + 2d)2(3 + 2^)^(5 + 2^)^(7 + 2d)2 and eigenvalues of T^^ will be ±iJ,k = ±Wk, with < < 1, 
B = d{A + d). We obtain ix^ by solving The eight for fc = 1, 2, 3, 4. ' 
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